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We study which and how a periodic orbit in phase space hnks to both the largest Lyapunov 
exponent and the expectation values of macroscopic variables in a Hamiltonian system with 
many degrees of freedom. The model which we use in this paper is the discrete nonlinear 
Schrodinger equation. Using a method based on the modulational estimate of a periodic 
orbit, we predict the largest Lyapunov exponent and the expectation value of a macroscopic 
variable. We show that (i) the predicted largest Lyapunov exponent generally depends on the 
periodic orbit which we employ, and (ii) the predicted expectation value of the macroscopic 
' variable does not depend on the periodic orbit at least in a high energy regime. In addition, 

. the physical meanings of these dependencies are considered. 

d 

§1. INTRODUCTION 



We are interested in characterizing chaotic Hamiltonian systems. One of the most basic 
indicators to characterize chaos for physicists is the largest Lyapunov exponent, and then 
Q> \ there have been attempts to estimate the largest Lyapunov exponents for several systems. For 

. low-dimensional Hamiltonian systems, Chirikov has found that the largest Lyapunov exponent 

, is close to an averaged eigenvalue in a strong chaotic regime for the standard map.^^ For 

I high-dimensional Hamiltonian systems, an analytical method to estimate the largest Lyapunov 

Q ■ exponents has been developed and applied. The method is based on a random approximation 

. and Riemannian geometry. For the Hamiltonian mean field model, "^^ which shows a second- 

order transition, the largest Lyapunov exponent has analytically been calculated. This study 
has shown that the relation between the largest Lyapunov exponent and the second-order phase 
transition. There are similar studies for the a-XY model, which is one of the extended models 



5— ( ' to study how non-additivity affects the statistics and dynamics. For the q-XY model, it has 



been shown that the largest Lyapunov exponent is a function of the interaction length. ^^^^^ 

It is important to estimate the expectation value of a macroscopic variable for a system 
with many degrees of freedom. Note here, the term "a macroscopic variables" refers to a 
quantity obtained by taking the average over many degrees of freedom. For instance, well- 
known macroscopic variables are the temperature and magnetization. In the equilibrium state, 
the expectation values of macroscopic variables can be estimated using tools of equilibrium 
statistical mechanics with some basic assumptions. These assumptions, such as ergodicity, the 
principle of equal weight, and so on, are not usually proved. Then estimating the expectation 
value of a macroscopic without such classical tools sheds light on basics of statistical mechanics. 

In dynamical systems theory, not only for Hamiltonian systems but for dissipative systems, 
utilizing and searching for periodic orbits are one of the most fundamental issues. When a 
periodic orbit is found in a system, we can see a part of the skeleton of phase space because the 
periodic orbit forms an invariant subset in the phase space. ( see Ref.,^-* for the Fermi-Pasta- 
Ulam-/3 model as an example ). 

There is the question whether a periodic orbit is related to the largest Lyapunov exponent 
and the expectation values of macroscopic variables. For the Navier-Stokes equation, it has 



recently been shown that values of some macroscopic quantities can be estimated using only 
one unstable periodic orbit. ^"^^^^^^ In order to clarify the meanings of these numerical studies, it 
has theoretically been studied that only one unstable periodic orbit is enough to derive relevant 
statistical property in a hyperbolic chaotic system with many degrees of freedom. In addition, 
a method to evaluate the correlation function has been proposed using both one periodic orbit 
and the projection operator method. Independently of such studies, it has been shown that 
a method with only one periodic orbit can predict the largest Lyapunov exponent in a class of 
models including the Fermi-Pasta-Ulam-/3 model. In Ref.,^^^ they have proposed the method 
predicting the largest Lyapunov exponent in a system with N degrees of freedom. The estimate 
consists of the following three steps: 

(1) Find one periodic orbit {q^'~' (t) , p^^ (t)) . 

(2) Estimate the linear growth rate along the periodic orbit, Tj {j = 1, ...,2iV), and define the 
instability entropy Sie ■= J2j(Tj>o) '^r 

(3) Predict Ai = 2Sie/N, based on the assumption Sje = Sks ~ AiA^/2, where Sks is the 
Kolmogorov- Sinai entropy. 

On the other hand, for a macroscopic variable 0{q,p), we could predict the expectation value 
in the equilibrium state. The method is given by the following procedure: 

• {0)po := 0{qPO{t),pPO(t))dt/To, where Tq is the period of 
This kind of substitution has also been used for the Fermi-Pasta-Ulam-/3 model in Ref.,^^^ where 
they have analytically shown that the largest Lyapunov exponent can also be obtained using 
an Riemannian geometric approach with this kind of substitution. 

If the above two estimates depend on the periodic orbit that we find, the dependency gives 
a clue to understand the link between averaged values (e.g., the largest Lyapunov exponent and 
the expectation value of a macroscopic variable) and a microscopic point of view (e.g, a periodic 
orbit). In general, it is difficult to discuss such a dependency, due to the lack of the number 
of periodic orbits. As a matter of fact, in the Fermi-Pasta-Ulam-/3 model, the number of such 
periodic orbits is five.-"^^^ To discuss the dependency, we need a chaotic model in which many 
periodic orbits are easily found. If the expressions of these periodic orbits are analytically 
written, these expressions give analytical expressions of linear growth rates along the orbits. It 
is noted here that there is a work in Ref.,^^) where they have compared with two periodic orbits 
in a different point of view. 

In this paper, we study the periodic orbit dependencies on (a) the method predicting the 
largest Lyapunov exponent^^) consisting of the three steps mentioned above, and (b) substitut- 
ing the expression of a periodic orbit into the definition of a macroscopic variable. We treat the 
discrete nonlinear Schrodinger equation, in which we exactly calculate the modulational esti- 
mates along N periodic orbits with N being the number of lattice sites. The fact that there are 
a number of analytically expressed periodic orbits is a feature peculiar to the discrete nonlinear 
Schrodinger equation. The existence of A*" exact periodic orbits allows us to discuss the periodic 
orbit dependencies for (a) and (b). Using these approach, we bridge some averaged quantities 
and a microscopic point of view: the largest Lyapunov exponent and the expectation value of 
a macroscopic variable from a periodic orbit. 

§2. Theoretical Prediction and Numerical Simulation 

The equations of motion and the Hamiltonian for the discrete nonlinear Schrodinger equa- 
tion are 



— = -t{uj+i + Uj-i - 2uj + -t\uj\ Uj) = — , (2-1) 
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with conditions n^+Tv = '^j, being the number of degrees of freedom. We restrict ourselves to 
the case that N is even. Here Uj are complex variables and 7 is a real parameter. The system 
(j2-ip has the conserved quantity / := X^jl^o^ l^jP' addition to the Hamiltonian. 

In the next subsection, we theoretically predict the largest Lyapunov exponent and the 
expectation value of a macroscopic variable. 

2.1. Theoretical Prediction 

First, we concentrate on the Lyapunov exponent, we predict the largest Lyapunov exponent 
with the three steps that are introduced in ^ 

As the step (1), we find the following periodic orbits in this model, 

^if°(t) = Afcexp{f(2^fci/iV-u;fc(|Afc|2)t)}, (2-2) 
^ki\Ak\^) ■■= -4.sm\7Tk/N)+j\Ak\\ 

where A; = 0,...,A — 1, and G C are the amplitudes of the orbits. It is noted that we do not 
make any approximation to express these unstable periodic orbits. 

As the step (2), we calculate the growth rate along u^^{t). To calculate the growth rates 
of the orbits, we substitute Uj = u^'^{l + //j), ^ 1 into Eq. ()2-ip . where Hj are complex 
variables describing the tangent phase space along u^^{t). After taking into account the linear 
terms only and using the Fourier transformation, we obtain the linearized equations around 
u^^{t) in Fourier space. 



Here 



dt V f^l-m J \ iB J\ f,l_^ 



N-1 N-1 
i — ^ /^r i 



(2-3) 



^i2TTmj /N 

j=0 ^ ^ ' m=0 

where m = 0, — 1, and B,C £M. are given by 

B = 4{sin2(7rA;/A) - sin2(7rm/A)} + j\Ak\ 



The growth rates of the periodic orbits (j2-2p are obtained as the eigenvalues of the linearized 
equation (|2-3p . The squared eigenvalues are calculated as follows, 

T(K:m] = — lb sm ( — ) — sm ( 

^ ' V N ^ ' 

-87|^.p(sin2(f )-sin2(^)). (2-4) 

Then we have the pairs of eigenvalues T±{k; m). It should be noted that the stability analysis has 
also been done without any approximation. The expression of eigenvalues shows the following 
symmetries for T±{k]m), 

T±{k;k) = 0, 

T±{k, N/2 -m) = T±{k, N/2 + m), 

T±{N/2-k,m) = T±iN/2 + k,m). (2-5) 
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When T^(A;; m) is positive, the periodic orbit labeled by k is unstable. On the other hand, when 
T'^{k;m) is negative, the periodic orbit is stable. Substituting Eq. (j2-2p into Eq. (j2-ip . we have 
the relation between the absolute value of the amplitude and the energy density e := H/N, {7i = 

As the step (3), the linear stability analysis gives an analytical prediction of the largest 
Lyapunov exponent, 

2 ^'"^'^ 

Ai = - ^+ik;m), (2-6) 

m=0(T2(fc;m)>0) 

and this can be viewed as the average of the growth rates. Here this expression of Ai depends 
on both the parameter k and the energy density e. 

Second, we predict the time-averaged value of a macroscopic variable using a substitution 
of the periodic orbit. This substitution is introduced in 311 As one of the macroscopic variables, 
we take 

^ _ 1 ^1 duj du* 
j=0 

Substituting the periodic orbits into Eq. (|2-7p . we predict the time-averaged values of the macro- 
scopic variable S. Then we have {S)po = In the high energy limit, we predict 
{S)po ~ 7^|Afc|^ ~ (— 2e7) provided ej < 0, and {S)po does not depend on k in this limit. In 
the low energy limit, there is a periodic orbit dependency. 

2.2. Numerical Simulation and discussion 

Let us compare our predictions with numerical simulations. We restrict ourselves to the case 
7 < 0. Numerical integrations of the canonical equation of motion are performed using a second 
symplectic integrator. ^^^^^^^ The time step of the integrator is set at 0.005 and it suppresses 
the relative energy error maxo<t<5oooo[{(-ff (i) — -f^(0)}/-?^(0)] ~ 10~^. Our initial conditions 
of Uj{0) are as follows. For Re{uj{0)}, small perturbation terms are added to oc cos(7rj), and 
lm{uj{0)} are exactly zero. The amplitudes of Uj{0) determine the value of the energy. The 
computing time is set to more than 50000 so that the time-averages converge, which are to 
obtain the largest Lyapunov exponent and the expectation value of the macroscopic variable. 
We do not take any ensemble average to obtain numerical data. 

First, we study the largest Lyapunov exponent. In Fig. [H we compare the largest Lyapunov 
exponent Ai obtained both the predictions (j2-6p with the numerical calculations for the model 
with N = 512. This figure shows that the analytical estimate depends on k, and that the 
prediction by taking k = 116 is best fitted to the numerical data in the regime 0.1 < E/N < 100. 
We denote it by A;*, and the value of A;* gives k^/N ~ 0.227 for = 512. For A^ = 128, the 
best fitting parameter is A* = 28 (i.e., /c^/A^ ~ 0.219), as shown in FigHJ Although we do not 
show any figure for cases A^ < 128, we find the non-trivial rule k^/N ~ 0.22. In Fig. [21 the 7 
dependency is shown for models with A^ = 128. For all 7 which we study, A;^, are the same each 
other. Then the rule k^/N ~ 0.22 can be applied for a wide range parameter regime. 

Let us look for the origin of the non-trivial rule k^/N ~ 0.22. Here we attempt to bridge 
the rule for the choice of the periodic orbit and components of the Lyapunov vector. Due to 
the definition of the largest Lyapunov exponent, the origin of the rule k^/N ~ 0.22 should be 
discussed in both the tangent phase space and phase space, not only in the phase space. When 
the rule k-^/N ~ 0.22 is applied to the periodic orbit (j2-2p . we can obtain the approximate 
periodicity in space, tij+4,5 ~ Uj. We identify these values, 4 ~ 5, with the correlation length. 
On the other hand, it is worth to note here that the localization phenomenon of components 
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Fig. 1. For the discrete nonlinear Schrodinger equation with A'^ — 512 and 7 = —1. Comparison of the analytical 
estimate with numerical data for the largest Lyapunov exponents. Pluses denote the numerically obtained 
the largest Lyapunov exponents and lines the predicted ones. The prediction using the periodic orbit labeled 
by A; = 116 is the best fit to the numerical data. 



10^ 
10' 
10° 
10-' 



10" 



I l| 1 TT] 1 TT] 1 TT] 1 TT] 1 FT^ 1 n 



+ 1^-0.1 + 

y=-l X 
■p-0.1 

- - - 



fr 

' ' 'I I l_J I l_J I l_J I l_J I L_J I l_L 



10° 10' 10^ 10^ lO'' 10^ 10^ 10^ 



E/N 



Fig. 2. For the discrete nonlinear Schrodinger equation with 7 = —0.1, —1, —2 and A'^ = 128. Comparison of the 
analytical estimate with numerical data for the largest Lyapunov exponents. Pluses, crosses and asterisks 
denote numerically obtained the largest Lyapunov exponents and lines the predicted ones with fc = 28. 
Predictions by periodic orbits labeled by = 28 are best fit to the numerical data. 



of the Lyapunov vector has been studied in Ref.^^^ In this paper, we calculate the normalized 
correlation length of components of the Lyapunov vector numerically as shown in Fig. [3j Here 
we define the correlation function for components of the Lyapunov vector as 

CO') = lim /* d^\du^{i^)dul{^)f. (2-8) 

t-^oo t - to 

Here to is taken as 100 for our calculations, duj{t) denote the j-th component of the tangent 
vector associated with the orbit Uj (t) . Fitting the normalized correlation function to a sum of a 
Lorentzian curve and an exponential curve, we extract the correlation length from the tangent 
space of the phase space. The fitting function which we use here is then 

where E, F and J are the fitting parameters. The value of this function equals unity at j = for 
any values of E, F and J. Due to this fitting function, we have two typical length scales. One of 
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Fig. 3. For the discrete nonlinear Schrodinger equation with = 128 and 7 = —1. Fitting the normalized 
correlation functions for the Lyapunov vector ( see Eq. (|2-8p ) to a sum of a Lorentzian curve and an 
exponential curve. A fitting curve gives the two correlation lengths, one of them is from the exponential 
curve and the other is from the Lorentzian curve. The correlation length is estimated as 4.6 from the 
exponential function for the high energy regime. 

them is from the width of the Lorentzian curve, -T, and the other is from the exponential curve, 
J. The correlation lengths obtained by J are approximately 4.5 ~ 5 for the large energy density 
regime. For example, the correlation length is 4.6 for the model with N = 128 and 7 = — 1 in 
the high energy regime as shown in Fig. [3l We can say that the correlation length in the tangent 
dynamics can be estimated using J. For the system with N = 128 and 7 = —0.1, although the 
value of J is in between 4.5 ~ 5 in the high energy regime, such correlation length is far from 
it, about 12 in the low energy regime. Then this explanation for the rule kt^/N ~ 0.22 is valid 
only in the high energy regime. Combining the consideration of components of the Lyapunov 
vector and that of the approximate spatial periodicity by applying the rule k/N ~ 0.22 to the 
periodic orbit, we could say that the emergence of the non-trivial rule is from the localization 
of the tangent dynamics in the high energy regime. It could be one of the reasons why the rule 
k-^/N ~ 0.22 appears, at least in the high energy regime. 

Let us consider the validity of the definition of the instability entropy. If we define another 
version of the instability entropy, S'j^ = X]m=o (A2(fc m)>o) -^+(^'"'')' ^he largest Lyapunov ex- 
ponent could be predicted as A'^ = {2/N) Ylm=o{x'^{k-m)>o) ^^+(^5 "^)' making the assumption 
that the Kolmogorov Sinai entropy is equal to Sj^. This new expression A'^, which is different 
from Eq. (j2-6p . cannot predict numerical data for the largest Lyapunov exponents in the whole 
energy density regime (no figure given). When we look for the parameter k^ which gives the 
best fit to the numerically obtained data in the high energy regime, we find k'^/N ~ 0.16 for 
the models with N = 512, 128 and 7 = -1. 

Next, we compare the prediction for the time-averaged value of S with numerical simulation. 
Fig. m shows that, in the high enough energy regime, each prediction labeled by k is in good 
agreement with numerically obtained data. In contrast to the success of the prediction in the 
high energy regime, it is difficult to conclude what happens in the low energy regime. Although 
we do not show any figure, in the low energy regime E/N < 100, the time-averaged values of 
S do not convergent fast in time. However, from the edge of the low energy regime, E/N ~ 50 
in Fig. m even the best fitted prediction, labeled by A; = N/2, differs from the data obtained 
numerically. Similar tendency is observed for the cases N < 512. 

Let us look for the reason why every periodic orbit (|2-2p can predict the time-averaged 
values of S in the high energy regime. Because the interaction range of this system is short, 
the system has additivity. For such an additive Hamiltonian system in a fully developed chaotic 
regime, we can define some macroscopic subsystems whose sizes are macroscopically arbitrary. 
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Fig. 4. For the discrete nonlinear Schrodinger equation with A'' — 512 and 7 = —1. Comparison of the analytical 
estimate with numerical data for S ( defined in Eq. (|2-7|) '). Pluses denote the numerical data and the lines 
theoretical ones. In the high energy regime, the theoretical estimate by every periodic orbit can predict data 
obtained numerically. However, in the low energy regime, the periodic orbit labeled by fc = N/2 is the best 
fit to the numerical data, and is far from the numerically obtained data. 



This kind of arbitrariness implies the suppression of the typical length scale, and then we can 
say that the expectation value of S does not depend on the periodic orbit characterized by 
k. It should be noted here that there is no contradiction between the /c-dependency for the 
prediction of the largest Lyapunov exponent and of the expectation value of the macroscopic 
variable. This is because the largest Lyapunov exponent is discussed in both the tangent phase 
space and the phase space, on the other hand, the expectation value of the macroscopic value 
is done only in the phase space. 

We conclude that (i) the prediction of the largest Lyapunov exponent depends on the 
periodic orbit, and (ii) the time-averaged value of the macroscopic variable can be predicted by 
any unstable periodic orbit, at least in the high energy limit. 

§3. CONCLUSIONS 

We have studied which and how periodic orbits predict both (a) the largest Lyapunov 
exponent and (b) the time-averaged value of a macroscopic variable in a Hamiltonian lattice. 

To clarify these questions, we have studied the discrete nonlinear Schrodinger equation in 
which there are N analytically expressed periodic orbits. In the nonlinear Schrodinger equa- 
tion, we have exactly constructed the modulational estimates along the A'^ periodic orbits and 
compared with numerical simulation. Then we have observed that the analytically predicted 
largest Lyapunov exponent depends on the periodic orbit and that there is a suitable periodic 
orbit for the prediction. The reason has been discussed in the phase space and the tangent 
phase space by studying components of the Lyapunov vector. On the other hand, to predict the 
time-averaged value of a macroscopic variable, we have substituted the analytical expression of 
the orbits into the definition of the macroscopic variable. The prediction is in good agreement 
with numerically obtained data in the high energy limit. The reason has been discussed in the 
phase space by considering the range of interactions. 

In this paper, we have focused on a way to bridge between the dynamics and statistics. 
Further investigations are necessary to quantitatively understand why the special periodic orbit 
can only predict the largest Lyapunov exponent even in a relatively low energy regime, and every 
periodic orbit can be used for predicting the time-averaged values of macroscopic variables in 
the high energy regime. We believe that this kind of investigations can help to elucidate the 
study of Hamiltonian systems with many degrees of freedom. 
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